Localization transition of stiff directed lines in random media 
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We investigate the localization of stiff directed lines with bending energy by a short-range random 
potential. Using perturbative arguments, Flory arguments, and a replica calculation, we show that a 
stiff directed line in 1+d dimensions undergoes a localization transition with increasing disorder for 
d > 2/3. We demonstrate that this transition is accessible by numerical transfer matrix calculations 
in 1 + 1 dimensions and analyze the properties of the disorder-dominated phase. On the basis 
of the two-replica problem, we propose a relation between the localization of stiff directed lines 
in 1 + d dimensions and of directed lines under tension in 1 + 3d dimensions, which is strongly 
supported by identical free energy distributions. This shows that pair interactions in the replicated 
Hamiltonian determine the nature of directed line localization transitions with consequences for the 
critical behavior of the Kardar-Parisi-Zhang (KPZ) equation. Furthermore, we quantify how the 
persistence length of the stiff directed line is reduced by disorder. 
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Directed lines (DLs) in random media or the more gen- 
eral problem of random elastic manifolds are one of the 
most important model systems in the statistical physics 
of disordered systems 1,. DLs in random media are 
related to important non-equilibrium statistical physics 
problems such as stochastic growth, in particular the 
KPZ equation [2], Burgers turbulence, or the asymmetric 
simple exclusion model (ASEP) [3 . Furthermore, there 
are many applications such as kinetic roughening [5] , pin- 
ning of flux lines [HE], random magnets, or wetting fronts 
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DLs in D = 1 + d-dimensional random media exhibit 
a disorder-driven localization transition for dimensions 
d > 2, which has been studied numerically for dimensions 
up to d = 3 [BHH]- In the context of the KPZ equation, it 
is a long-standing open question whether there exists an 
upper critical dimension, where the critical behavior at 
the localization transition is modified (for a recent discus- 
sion, see, e.g., [HI])- We will address this question from a 
new perspective by revealing a relation between DLs in 
l + 3d dimensions and stiff directed lines (SDLs) in 1 + d 
dimensions in a random medium, which we validate for 
d = 1 by numerical simulations. We define SDLs as lines 
with preferred orientation and no overhangs with respect 
to this direction, but a different elastic energy penaliz- 
ing curvature instead of stretching, cf. Fig. [T] We inves- 
tigate their disorder-induced localization transition and 
the scaling properties of conformations in the disordered 
phase. Because the proposed relation between DLs and 
SDLs is based on return probabilities of pairs of replicas 
our results suggest that the critical properties of DLs in 
a random potential and, thus, of the KPZ equation are 
determined by the corresponding two replica problem. 

Model of the stiff directed line. The configuration of 
a directed line, i.e. one without overhangs or loops and 
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FIG. 1: (color online) Paths with lowest energy for given end- 
ing states (top, globally optimal path thicker) for (a) the stiff 
and (b) the tense directed line (right). The dots mark favor- 
able regions of the Gaussian random potential V (realizations 
of the quenched disorder are not identical in a and b). 



without inextensibility constraint, can be parametrized 
by (x,z(x)) (0 < x < L) with a d-dimensional displace- 
ment z(x) normal to its preferred orientation. Each con- 
figuration of a SDL is associated with an energy 



H 



dx 



'-(d 2 x z(x)) 2 + V(x,z(x)) 



(1) 



where the first term is the bending energy (to leading 
order in z). The second term is the disorder energy 
with a Gaussian distributed quenched random potential 
V(x,z(x)) with zero mean V — and short-range corre- 
= g 2 5 d (z - z')S(x 



lations V(x, z)V(x', z') = g 2 S d (z - z')5(x - x 1 ) [35]. 

The SDL model ([I]) is an approximation to the so- 
called worm-like chain or Kratky-Porod model [HI [12] 
Uwlc = Jo ds [§(<9j?r(s)) 2 + V~(r(s))] which is a com- 
mon model for semiflexible polymers, such as DNA 
or cytoskeletal filaments like F-actin. The chain is 
parametrized in arc length, leading to a local inexten- 
sibility constraint |<9 s r(s)| = 1, which is hard to account 
for, both numerically and analytically jT3l [14] . The ap- 
proximation ([I]) only applies to a weakly bent semiflexible 
polymer on length scales below the so-called persistence 
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length, which is 36J L p = (D — 1)k/2T for thermally fluc- 
tuating semiflexible polymers [14] . Above the persistence 
length, a semiflexible polymer loses orientation correla- 
tions and starts to develop overhangs. Also in a quenched 
random potential the SDL model describes semiflexible 
polymers in heterogeneous media, only as long as tan- 
gent fluctuations are small such that overhangs can be 
neglected, which is the case below a disorder-induced per- 
sistence length, which we will derive below. 

We consider the SDL model also in the thermodynamic 
limit beyond this persistence length, because we find ev- 
idence for a relation to the important problem of DLs 
in a random medium in lower dimensions. This relation 
is based on replica pair interactions and shows that pair 
interactions also determine the nature of the DL localiza- 
tion transition. Moreover, this relation can make the DL 
transition in high dimensions computationally accessible. 
We will now outline the idea behind this relation. 

Relation to directed lines. The difference between 
SDLs and DLs [I] is the second derivative in the first 
bending energy term in eq. ([!]) for SDLs, which differs 
from the tension or stretching energy ~ J dx(d x z(x)) 2 
of DLs. This results in different types of energetically 
favorable configurations (see Fig. [TJ: large perpendicu- 
lar displacements z are not unfavorable as long as their 
"direction" does not change. Displacements z are charac- 
terized by the roughness exponent £, which is defined by 
(z 2 (L)) ~ L 2C -. The thermal roughness is Cth,r = 1/2 for 
DLs and (th.n = 3/2 for SDLs. Here and in the following 
we use subscripts r (tension) and k (bending stiffness) to 
distinguish between the two systems. 

Although typical configurations are quite different, a 
SDL subject to a short-ranged (around z = 0) attrac- 
tive potential V(z) can be mapped onto a DL in high 
dimensions d' — 3d [TS1 US]. This equivalence is based 
on the probability that a free line of length L starting at 
(z(0) = 0) "returns" to the origin, i.e., ends at z(L) = 0. 
This return probability is characterized by a return ex- 
ponent x- Prob(z(L) = 0) ~ L~ x . For DLs, which are 
essentially random walks in d transverse dimensions, the 
return exponent is \ T = d/2 [17] . whereas it is \k = 3d/2 
for a SDL (after integrating over all orientations of the 
end) [TS] ; they are related to the roughness exponents by 
X = (,d [HI- The return exponent governs the binding 
to a short-range attractive potential or, equivalently, the 
binding of two lines interacting by such a potential |15l - 
117] . The relation Xr(3d) = Xn{d) implies that the bind- 
ing transition of two DLs in 3d dimensions maps onto to 
the binding transition of two SDLs in d dimensions. 

In the replica formulation of line localization problems 
such as ([!]), the random potential gives rise to a short- 
range attractive pair interaction (see below). Further- 
more, the critical temperature T C T for a DL in a random 
potential is believed to be identical to the critical tem- 
perature T2 jT for a system with two replicas [9] [19]. One 
aim of this work is to generalize this conjecture by demon- 
strating that the d —> 3d analogy between DLs and SDLs 
in a random potential holds for the entire free energy dis- 
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FIG. 2: Sketch of the phase diagram as implied by Flory- 
type arguments. The arrows indicate the flow of the disorder 
"strength" g under renormalization. 



tribution. Because this analogy is rooted in the binding 
transition of replica pairs, we can conclude that critical 
properties of the localization transition are determined 
by pair interactions in the replicated Hamiltonian, which 
has been previously suggested in Refs. [20112I]- Moreover, 
it has been proposed that pair interactions can be used 
to formulate an order parameter of the disorder-driven 
localization transition of DLs in terms of the overlap 

q = Iimi^oo lf£dx8(zi(x) - z 2 (x)) [22H2H, Le -> the 
average number of sites per length, that two lines in the 
same realization of the disorder have in common. Local- 
ization by disorder gives rise to a finite value of the pair 
overlap q. 

Scaling analysis. We start with a scaling analysis of 
SDL behavior in disorder by a Flory-type argument. For 
displacements ~ z the bending energy in eq. (nj) scales 
as Ef, ~ kz 2 L~ 3 , which also leads to (z 2 ) ~ Lr/L v and 
the thermal roughness exponent Qth.n = 3/2. The dis- 
order energy in eq. |lj scales as E4 ~ g\/ Lz~ d . Using 
the unperturbed thermal roughness in the disorder en- 
ergy we get E d ~ L( 2 " 3d )/ 4 , from which we conclude 
that the disorder is relevant below a critical dimension 
d < rf c ,K = 2/3. For d > d CyK — 2/3 and, thus, in all 
physically accessible integer dimensions, the SDL should 
exhibit a transition from a thermal phase for low g to a 
disorder dominated phase above a critical value g c of the 
disorder (see Fig. [2]). In the disorder dominated phase, 
the SDL becomes localized and assumes a roughened con- 
figuration, see Fig. [I] 

Balancing the Flory estimates, Ef, ~ Ed, gives the 
roughness z ~ L^ Fl , if disorder is relevant. This leads to 
£fi,k — 7/(4 + d), which is applicable below the critical 
dimension d < d CiKl where (fi,k > Cth,K- Above the crit- 
ical dimension, the Flory result would give Cfi,k < (,th,K, 
which contradicts a roughening of the SDL as it adjusts 
to the random potential. Furthermore, the exponent u> 
related to the energy fluctuations via A_F ~ would be 
negative, since the scaling of the energy implies a general 
scaling relation u) = 2£ re — 3 (note that we do not subscript 
w as we believe w T = ui K , see below) . An exponent uj < 
contradicts the existence of large disorder-induced free 
energy fluctuations in the low-temperature phase [25]l26j. 
for which there is also strong numerical evidence [7H5]. 
We conclude that this kind of argument is not applicable 
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above the critical dimension. The same problems occur 
in Flory arguments for DLs for d > d CtT = 2 as well as in, 
for example, functional renormalization group analysis 
[27] for d 2.5. We conclude that, in contrast to DLs, it 
is sufficient to study 1+1 dimensional SDLs in numerical 
transfer matrix calculations in order to explore the prop- 
erties of their localization transition such that we focus 
on this case in the following. 

Variation in replica space. To go beyond scaling 
arguments we use the replica technique |28j follow- 
ing the treatment of directed manifolds [22], the re- 
sults of which we summarize briefly in the follow- 
ing. We write the averaged and replicated parti- 
tion function as Z n — Yl a (JVzoc)exp(—j3H T ep) with 
the replica Hamiltonian in Fourier space "Hrop = 

^YL^ki^+riZ-^TZ^ffah^) where 
- z/3 and with an attractive potential f\(z) 
of range A. As mentioned before, T-L^ is related to a 
pair binding problem: in the limit A ~ the second term 

becomes J2 a ,p Jo dx5(z a - zp). 

Using variation in replica space we find one-step replica 
symmetry breaking for d > c? CjK and can show that there 
is no localized solution unless the potential strength g and 
correlation length A are above finite values. We interpret 
this as an indication for the existence of a critical disorder 
strength or a critical temperature for d > d C K . 

Numerical results. More progress is possible by exten- 
sive numerical studies using the transfer matrix method 
[H |3U] both for T = (see Fig. ^ and for T > 0. The 
transfer matrix element for a segment of a SDL with 
length AL = 1 starting at z with orientation dz/dx = v 
and ending at z' with v' follows from a discretization of 
(flj) as AE — k/2(v — v') 2 + V(z'), where we have chosen 
z = z + v' to simplify the transfer matrix. For the sake 
of simplicity we choose g = 1 and vary the temperature. 

In order to test the relation between DLs and SDLs, 
we compare free energy fluctuations for DLs in 1+3 and 
SDLs in 1+1 dimensions. Previous direct numerical stud- 
ies of DLs in 1+3 dimensions found an exponent w « 0.18 
[HI [5] in the low temperature phase, the most precise value 
from kinetic roughening studies is lu s» 0.186 [31]. We 
compare these values with own numerical transfer matrix 
calculations for SDLs in 1+1 dimensions. We determine 
the exponent lu directly by fitting AF — (F 2 — F 2 ) 1 / 2 cx 
L u , which gives values w « 0.16 for DLs in 1+3 dimen- 
sions, or by studying the distribution of the free energy 
shown in Fig. [3ji, which is obtained by computing the free 
energy for every sample and rescaling to zero mean and 
unit variance, G F {x) = Prob((F - F)/AF = x). The 
asymptotic behavior of the negative tail of the rescaled 
free energy distribution for low temperatures, which is 
of the form lnG F (x) — -\x\ ri (x < 0,\x\ > 1), allows 
us to determine the energy fluctuation exponent lu via 
the Zhang argument pQ giving r] = 1/(1 — w). We get 
rj 1.23 and therefore u) f=s 0.19. These values agree 
with the literature values for DLs in 1+3 dimensions. 

For direct comparison of the entire free energy distri- 



butions of a SDL in 1+1 and a DL in 1+3 dimensions 
we simulated both systems and find that the rescaled free 
energy distributions in the low temperature phases have 
to be considered identical within numerical accuracy, see 
Fig. [3^,. This strongly supports the relation between DLs 
in 1 + 3d dimensions and SDLs in 1 + d dimensions. 

In a third approach, we can calculate £ K and lu = 
2£ K — 3 by measuring a "local" version of the rough- 
ness exponent [32J 2((L) = \og 5 {z 2 {L)/z 2 (L/5)). The 
data shown in Fig. [3Jd shows two distinct regimes for 
high and low temperatures and a significant "dip" around 
T = 1.4. For low temperatures, values 2£ K w 3.11 are 
consistent with uj ps 0.11. As for DLs [33], it can be ar- 
gued that to should vanish at the transition, resulting in 
fluctuations of the free energy that scale logarithmically 
with L, AF ~ In 1 / 2 L [9], resulting in a roughness ex- 
ponent Ck = 3/2. This seems to hold, even though the 
numerical value for high temperatures is slightly above 
f re = 3/2, which is strong evidence for a phase transition 
at T c w 1.4. We support this by studying the difference 
of quenched and annealed free energies SF = F — F ann , 

1/2 

which should as well scale as 8F ~ In ' L at the transi- 
tion. We determined -Fann by simulating a system with- 
out disorder and adding the contribution of the annealed 
potential, F ann = F g=Q - Lg 2 f3/2, see Fig. ^p. 

Finally, we identify an order parameter of the localiza- 
tion transition. For DLs, the disorder-averaged overlap 

q = linif^oo jr dx5{z\(x) — z 2 (x)) of two replicas has 
been proposed as order parameter [22J 23J . Up to now, it 
has been numerically impossible to verify this order pa- 
rameter for DLs in d > 2 dimensions where a localization 
transition exists because the relevant 2<i-dimensional two 
replica phase space is too large. For SDLs, on the other 
hand, the transition is numerically accessible already in 
1+1 dimensions and we show that the overlap q is indeed 
a valid order parameter using an adaptation of the trans- 
fer matrix technique from Ref. [53] > see Fig. This in- 
volves simulating two interacting SDLs, therefore we can 
only use lengths up to L = 30 and 10 3 samples. For DLs, 
it has been found that the overlap at criticality decays 
as q - L s with S = -2( = -(1 + oS) in d = 3 [52] . 
This has been extended to finite temperatures yielding 
q ~ \T — T c \~ uTl . Indeed, we find a qualitatively similar 
behavior q ~ |T — T c |~^ with an exponent f3' s» 1.3 — 1.4. 
Our best estimate for S is E ps —0.75. For the correla- 
tion length exponent v we find values v ~ 2 compatible 
with the corresponding problem of DLs[HJ [U [T^; such 
that our present results deviate from j3' = vll. Because 
of small simulation lengths L we do not conclude this 
deviation to be a definite statement against the renor- 
malization group results presented in Ref. [52]. Never- 
theless, the connection between DLs and SDLs provides 
the first system to test the proposed order parameter in a 
localization transition numerically and to determine the 
otherwise inaccessible exponents /3' or S. 

Disorder-induced persistence length. As stated be- 
fore, the model of a SDL in a random potential de- 
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FIG. 3: (color online) a): Rescaled G F {X) = P((F - F)/AF) free energy distribution for a stiff directed line (SDL) in 1 + 1 
dimensions. We show distributions for T — (light (blue) squares, ground state energy) as well as for three finite temperatures 
T < T c (black thin solid line), T w T c (red, thicker solid line) and T > T c (green, light thin solid line). Results for a 
directed line (DL) in 1+3 dimensions are shown (dark circles), b): Local (see text) roughness exponent 2£ K as a function of 
T. Deviations from the analytical value 2( K = 3 at high temperatures indicate numerical problems, nonetheless there is a clear 
"dip" at T w 1.4, which we identify as the critical temperature. For low temperatures, we find values consistent with cj w 0.11. 
Inset: Reduced free energy SF = F — F ann rescaled by In 1 / 2 L for lengths L = 50, 60, . . . , 100 as a function of T, There is a 
pseudocrossing around T ps 1.38. c): The overlap order parameter q, as a function of T. We estimated q from finite lengths 
using a fit qr(L) = a(T)/L + g 00 (T). Inset: Double-logarithmic plot of the overlap q versus T c — T (with T c — 1.44), the solid 
line is given by q ~ (T c ~ T)~^' with w -1.36. 
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FIG. 4: Disorder-reduced persistence length L p of the SDL 
for g = 1. L p matches its thermal value (solid line) for T > T c 
and is reduced and approximately constant for T < T c . Inset: 
L p at T = versus the potential strength g. The Flory-result 
Lp ~ g" 1 (dotted line) matches the data. 



scribes a worm-like chain in a heterogeneous environ- 
ment on small length scales L < L p such that over- 
hangs can be neglected, which gives the defining cri- 
terion (v 2 )(L = L p ) = 1 for the persistence length. 
The line roughens in the low temperature phase, which 
gives rise to a persistence length decrease as compared 
to the thermal persistence length L p ~ k/UbT. At low 
temperatures, the Flory-result z ~ (g/ K) 2 ^ i+d ^ L 7 ^ 4+d ^ 
leads to L p ~ («/ 'g) 2 /^ d ) . For d = 1, we deter- 
mine L p numerically via the above defining criterion, and 
the data presented in Fig. [4] does indeed show an only 
weakly temperature-dependent disorder-reduced persis- 
tence length for T < T c . The T — results match the 
Flory-result L p ~ re/g. 

Conclusion. We studied stiff directed lines (SDLs) in 
1 + d dimensions subject to quenched short-range ran- 
dom potential analytically and numerically. Using Flory- 
type scaling arguments and a replica calculation we show 
that, in dimensions d > 2/3, a localization transition ex- 



ists from a high temperature phase, where the system is 
essentially annealed, to a disorder-dominated low tem- 
perature phase. The low temperature phase is charac- 
terized by large free energy fluctuations with an expo- 
nent to > 0, which cannot be calculated by scaling ar- 
guments or replica calculations. By extensive numerical 
transfer matrix calculations in 1+1 dimensions we find a 
value oj w 0.18, which is close to the established value 
uj fa 0.186 for directed lines (DLs) under tension in 1+3 
dimensions. Moreover, the rescaled free energy distribu- 
tions are identical. Both points suggest that the nature 
of the low-temperature phase is very similar, if not iden- 
tical. 

This strongly supports a relation between DLs in l+3c£ 
and SDLs in 1 + d dimensions, which is based on identical 
return exponents x f° r two replicas to meet. The validity 
of a relation based on properties of a single replica pair 
suggests that the critical properties of DLs in a short- 
range random potential are governed by replica pair in- 
teractions. The mapping can make DL transitions in 
high dimensions computationally accessible, which we 
demonstrated in showing that the two-replica overlap 
provides a valid order parameter across the localization 
transition of SDLs in 1+1 dimensions. Furthermore, the 
importance of pair interactions suggests that the criti- 
cal temperature for DLs in random potentials is indeed 
identical to the temperature below which the ratio of the 
second moment of the partition function and the square 
of its first moment diverges. The binding transition of 
DL pairs becomes discontinuous for d > 4 and, analo- 
gously, the binding of SDL pairs for d > 4/3 [TU [33]. 
Because DLs in random potentials are equivalent to the 
KPZ equation [2], the validated relation to the SDL sug- 
gests that the roughening transition of the KPZ problem 
could acquire similar discontinuous features for d > 4 
dimensions. Finally, we calculated the reduction of the 
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